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ABSTRACT 

Purely gravitational perturbations are considered in a thin rotating disk composed 
of several gas and stellar components. The dispersion relation for the axisymmetric 
density waves propagating through the disk is found and the criterion for the local 
axisymmetric stability of the whole system is formulated. In the appropriate limit 
of two-component gas we confirm the findings of Jog & Solomon (1984) and extend 
consideration to the case when one component is collisionless. Gravitational stability 
of the Galactic disk in the Solar neighborhood based on the multicomponent insta- 
bility condition is explored using recent measurements of the stellar composition and 
kinematics in the local Galactic disk obtained by Hipparcos satellite. 

Key words: Galaxy: evolution - Galaxy: kinematics and dynamics - Galaxy: solar 
neighborhood - ISM: kinematics and dynamics - stars: kinematics 



1 INTRODUCTION 

It was first shown by Safronov (1960) that a thin gaseous ro- 
tating disk can become unstable against local axisymmetric 
perturbations under the action of its own gravity. Quan- 
titatively the instability criterion is expressed in terms of 
Toomre's Q: 

(1) 
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with k being the epicycle frequency, c being the speed of 
sound, and E being the surface density of the disk. Instabil- 
ity arises when Q < 1. 

Disks composed of stars behave very similarly to the 
gaseous disks for long wavelength because epicyclic motion 
of stars is the same as that of the gas. But collisionless disks 
of stars have no pressure to resist short wavelength density 
perturbations and thus the density waves in this case are 
different from simple sound waves present in the gaseous 
disks in this regime (Toomre 1964). When the wavelength 
becomes shorter than the typical amplitude of the star's 
epicyclic motion the gravitational effect of the perturba- 
tion averages out and stars perform almost pure epicyclic 
motion with frequency k. Instability can occur though for 
some intermediate wavelength as in the case of the gaseous 
disks. Toomre (1964) had shown that in the stellar disks 
self-gravity can support exponentially growing axisymmet- 
ric modes when the condition Q < 1 is satisfied with n 
replaced by 3.36 and c replaced by the stellar velocity dis- 
persion a in the radial direction in definition (|l|) . Except for 
the qualitative difference in behavior at small wavelength, 
stellar and gaseous disks behave almost identically for large 



wavelength which leads to a close similarity of the corre- 
sponding instability conditions. 

However, real disks of spiral galaxies do not consist of 
only one component which could be characterized by a single 
value of the velocity dispersion. Indeed, a significant part 
of the ISM in galactic disks is usually concentrated in the 
form of cold gas with characteristic sound velocity ~ 5 — 10 
km s~ . At the same time stars usually have larger velocity 
dispersions ~ 20 — 35 km s - . 

The problem of the treatment of such two-component 
disks was first addressed by Jog & Solomon (1984a, JS). 
They assumed galactic disk to consist of two fluid compo- 
nents: one hot, playing the role of the stellar population, and 
one cold representing the usual ISM, and they neglected the 
fact that in real galaxies hot component is collisionless and 
provides no damping. Though this could be important short- 
coming of the analysis for the short wavelength it turns out 
to be a good approximation for the stability purposes as we 
will see later in this paper. Jog & Solomon (JS) have demon- 
strated that even relatively small amounts of the cold gas 
(~ 10% by mass) can very effectively destabilize the whole 
system. 

The two-fluid stability criterion was applied to a variety 
of disks by Jog & Solomon (1984b), Elmegreen (1995), and 
Jog (1996). They showed that our Galaxy is stable against 
local axisymmetric perturbations to the extent of the ac- 
curacy of the parameters and mass distribution models as- 
sumed. Theory can be further refined by including the effects 
of the finite thickness of the disk on its stability (Toomre 
1964) and by taking into account the vertical motions in the 
disk (Romeo 1992). It was shown that the instability con- 
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dition remains the same provided the unperturbed surface 
density is multiplied by an appropriate reduction factor. 

Real galaxies always contain more than two isothermal 
components. It is now widely believed that ISM in the galac- 
tic disk is subdivided into a number of components with 
different temperatures and, thus, different dynamical prop- 
erties (McKee & Ostriker 1977). The same is true about 
the stellar population. Stars in the galactic disk are being 
constantly scattered by giant molecular clouds and transient 
spiral arms which steadily increase their velocity dispersions 
(Spitzer & Schwarzschild 1951, 1953; Barbanis & Woltjer 
1967; Carlberg & Sellwood 1985). It means that the stars 
of different ages have different velocity dispersions and thus 
dynamically should be treated separately. 

Recent results of the Hipparcos satellite (ESA 1997) 
provided us with a wealth of information about the local 
stellar kinematics. The proper motions of different star pop- 
ulations were accurately measured and the resulting velocity 
ellipsoids were constructed for them enabling one to study 
different kinematic constituents of our Galaxy separately 
(Dehnen & Binney 1998; Mignard 2000). This mission also 
gave us a lot of information about distances to the stars 
which provided us with the local densities of various types 
of the stars (Holmberg & Flynn 2000). 

In this paper we present the analysis of the axisym- 
metric gravitational stability of the thin rotating disk which 
consists of a number of components, each of them being 
characterized by its own temperature and surface density. 
In doing so we distinguish between the two types of com- 
ponents: collisional and thus having pressure forces such as 
usual gas, or collisionless such as stellar component, which 
needs a kinetic treatment. The derivation of the dispersion 
relation is presented in §^ and|^. In §^ we compare our stabil- 
ity condition with the one derived in JS to test the difference 
arising when one of the disk components is treated as colli- 
sionless which better represents real galaxies. Finally, in §|| 
we apply our results to the case of our Galaxy in the Solar 
neighborhood using recent data from Hipparcos satellite. 



2 BASIC EQUATIONS 

We work in non-rotating cylindrical system, r, (j>, z, such that 
z-axis coincides with the rotation axis of the disk and angle 
<j> increases in the direction of rotation. We start first with 
the equations for the fluid components and then make some 
correction to take collisionless components into account. We 
will usually use index i to describe gas components and index 
j for stellar components. 

All the components are spatially concentrated in the 
thin disk and the effects of the finite disk thickness will usu- 
ally be disregarded. All the motions are assumed to occur 
only in the plane of the disk. 

We suppose that there are n g gaseous components con- 
tributing to the mass of the Galaxy, each characterized by 
the sound velocity d and surface density E 9 i, and n s col- 
lisionless components with E SJ being the surface density of 
the j'-th component. Each collisionless component in the un- 
perturbed state is assumed to have a Schwarzschild distri- 
bution function, that is for j-th component 
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where v c is the circular velocity at current distance from the 
galactic center and a r j and a^j are the velocity dispersions 
in r and <f> directions correspondingly. These velocity dis- 
persions are related by <7 2 /a 2 = 4B 2 /k 2 , where B is Oort's 
B constant (Binney & Tremaine 1987). We will later use 
simply a to denote the velocity dispersion in the r-direction 

For the description of collisional components we assume 
usual hydrodynamical description with isotropic pressure. 
To describe the star-like components we use a different ap- 
proach which is rooted in the kinetic treatment of the colli- 
sionless systems as described elsewhere (Toomre 1964; Bin- 
ney & Tremaine 1987). We will see in §^ that this accu- 
rate treatment shows that fluid approach really does quite 
a good job in describing the stability of the collisionless sys- 
tems with multiple components as it does in one-fluid case, 
though some quantitative differences exist. 

Equations governing the motion of the gas components 
in our coordinates are Euler's and continuity equations: 



^i + V(E iVi ) = 0, 



(3) 
(4) 



for z-th component. 

We assume that pressure of each component Pi = KT. 
and introduce specific enthalpy 
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We linearize equations (^|) and (0) by assuming that 
Vri — Ui,v<f>i = v c + Vi,hi = hoi + hu, Si = Eoi + Eij, and 
$ = $o + 3>i- 

Then equations ^ and (Q) reduce in the first order to 
(Binney & Tremaine 1987) 



dt + 50 

dm dm 
dt d(j> 
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where 



B = 



d(Qr) 
dr 



(6) 
(7) 
(8) 

(9) 



is the Oort's B constant. 

Since here we are interested in axisymmetric perturba- 
tions only it is always supposed that d/d<f> = in our analy- 
sis. We assume that all the first order dependent quantities 
vary like exp[i(fcr — wt)], where ui is the angular frequency 
and k — 2ti/\ is the wavenumber. For the local analysis 
we employ the WKB approximation (or tight-winding ap- 
proximation) which requires that kr 1 and allows us to 
neglect terms proportional to 1/r compared to the terms 
proportional to k. With all these simplifications equations 
(H)-(H) reduce to 



2FlVi = — + hu), 
—iuiVi — 2Bui — 0, 
— iwEii + ikEoiUi = 0, 



(10) 

(11) 
(12) 



for i-th gas component. 
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These equations have to be supplemented with Poisson 
equation 



4nG 



(13) 



3=1 



where S(z) is the Dirac delta function arising from our as- 
sumption of infinitely thin disk. With and $1 being 
proportional to exp[i(fcr — ut)] the relation between the per- 
turbations of the potential and surface densities becomes 
(Toomre 1964) 




X>" + X> 



(14) 
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The sound speed for each gas component is defined as 
Cj — dPj/dSoj , and in this case the perturbation of enthalpy 
reduces to 



(15) 



3 DISPERSION RELATION AND STABILITY 
CRITERION 

Equations @, and @ form a closed set of lin- 

ear equations and we now solve them to get the dispersion 
relation. 

From equations (|l(j) and ( |i~ij ) we can relate perturba- 
tion of the radial velocity of each component Ui to the po- 
tential perturbation $i using equation (|Te|): 



.0± U 1+ c 2 ^ 



where 



(16) 



(17) 



and k 2 = — AQ.B is the epicyclic frequency. 

Now we are in position to treat the star-like compo- 
nents. Since they are supposed to be collisionless they cannot 
create any pressure and it means that for j'-th stellar com- 
ponent sound velocity Cj — 0. But because of the epicyclic 
motion at any given point in the disk there are stars from 
different parts of the perturbed structure and it leads to im- 
portant cancellation effects. They were first calculated by 
Toomre (1964) and it can be shown that each stellar com- 
ponent is described by the Jeans equations which are pretty 
similar to the usual hydrodynamic equations with important 
difference being that instead of equation (n6|) we have 



cok . _ / lo k a a 
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where the reduction factor F is given by the expression (Bin- 
ney & Tremaine 1987) 



F{s,X) = -(1 - « )e" 
X 



Ux) 

— s 2 /n 



(19) 



and I n are the Bessel functions of order n. 

Now, from @, © and © we can eliminate m to 

get 



k 2 E 0i 



A + k 2 c 2 
for i-th gas component and 

7 2n 



(20) 



(21) 



for j-th stellar component and Tj = J-(lo/k, k 2 a 2 /k 2 ). 

We can then substitute ( |2"l| ) and (|2^) into ( pr| ) to obtain 
finally the desired dispersion relation: 
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E 
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(22) 



If only 2 fluid components are present in the disk this 
dispersion relation clearly reduces to the one derived by JS. 
If one of the components is collisionless the dispersion rela- 
tion is different and it was first considered by Romeo (1992) . 
We discuss the difference between the two cases in «. 

The multicomponent disk is unstable if lo 2 < because 
in this case oscillatory behavior of the density waves changes 
to an exponential growth. The dispersion relation (^) has 
an infinite number of solutions with u> 2 (k) > for all k, 
which are clearly stable. But there is also one mode of os- 
cillations in the system which has a single solution with 
— oo < lu 2 I k 2 < 1 and this mode could be unstable. 

To find the condition for instability we note that as 
uj 2 — » — oo at a fixed wavenumber k LHS of the equation 
( p2| ) tends to zero, and it monotonically increases to +oo 
as uj 2 — * k, 2 . This property of a steady growth is important 
because if ( p2|) has no solution for negative lj 2 it implies that 
the LHS of (p3) at uj 2 = must be less than 1 and vice versa. 
So the axisymmetric instability arises in multicomponent 
disk if and only if 



3=1 



where 



2 / k 2 a 2 
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Since 

n=l 

(Abramowitz & Stegun 1964) we can rewrite (S) as 



9 4 



1 - e 



1 I 



I (k 2 a 2 /K 2 ) 



k 2 a 2 j^ 2 



(23) 



(24) 



(25) 



(26) 



Suppose again that there are only two fluid (collisional) 
components in the disk then the instability condition ( |23[ ) 
reduces to the one derived in JS as it should. 



4 REALISTIC TWO-FLUID CASE 

Now we consider the stability in the case when we have only 
two fluids present in the disk in more details. Our main pur- 
pose here is to see how the fact that one of the components 
is in reality collisionless changes the overall dynamics of the 
system from the case considered by JS. 
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In the case of two components, one stellar and one 
gaseous, we assume E s and E s to be the unperturbed gas 
and stellar surface densities, c 9 the gas sound speed, and a a 
the stellar velocity dispersion in r-direction and introduce 
the following dimensionless quantities 



Q s = 



TvGT, a 
q = k<Jsj K, 



7rGE 



R = Cg/tTs 



(27) 
(28) 



Here Q g is a Toomre's Q parameter for gas, while Q s is 
different from the Toomre's Q parameter for the collisionless 
system, which is given by k<7 s /(3.36GE s ) (Toomre 1964). It 
means that with our definition of Q B the one-component 
stellar disk is gravitationally unstable when Q s < 3.36/-7T = 
1.07. 

With these definitions our instability condition (??) re- 
duces to 



2 

Q~s 



-R- 



> 1. 



(29) 



1 + q 2 R 2 

If one follows the Jog & Solomon approach and treats 
stellar component as a fluid with sound speed equal to a s , 
one gets the following instability condition in terms of our 
dimensionless variables 
2 



> 1. 



(30) 



Qsl + q 2 ' Q a ' l + q 2 R 2 

In Figure [l] we show the dependences of the uj 2 upon 
k for some particular choices of parameters and compare 
the curves in the case of gas-gas mix considered by Jog 
& Solomon (labeled as G+G) and star-gas mix (labeled as 
S+G), calculated using dispersion relation ( p2j ) for the case 
of two-component disk. If one adopts k = 36 km s _1 kpc -1 
to characterize the disk rotation and sound speed of the gas 
c = 5 km s _1 , then the model corresponding to the Figure 
[l|a has E s = 45 M pc -2 , E 9 = 9 M Q pc" 2 , and a = 25 
km s~\ which is a good representation of the typical spiral 
galaxy like our own. 

Second plate of the same Figure correspond to the 
model with E s = 17 M pc -2 , E s = 11 M pc -2 , and 
a — 12.5 km s _1 for the same choice of k and c. This might 
be thought of as a central part of a relatively young spiral 
and it turns out to be unstable with both instability criteria. 

One can immediately see that in general the difference 
between the curves given by conditions ( pS(i| ) and ( p9[ ) is small 
for low enough wavenumbers as it should be because pres- 
sure forces are negligible in this regime and it is only the 
epicyclic motions and the self-gravity which determine the 
dynamics of the disk. However, as k grows, in gas-gas case 
the frequency describing the oscillations of the whole sys- 
tem grows due to the pressure forces in the gas, while in 
the star-gas case this does not happen. Despite the presence 
of the gas in the system it is really stars which determine 
the natural frequency of oscillations of the system in this 
case and as k — > oo frequency becomes constant - lu — > k 
- and it presents a dramatic difference in the disk response 
to the axisymmetric density perturbations at small wave- 
lengths compared to the fluid-fluid case. 

In the two-fluid case JS found that the ratio of the per- 
turbation of the gas surface density to that of the stellar 
surface density is always a growing function of wavenumber 
k and for large k this ratio becomes very large. This is not 
the case when one of the disk components is collisionless. 
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Figure 1. Plots showing the dependence of the normalized square 
of the frequency u) 2 / k 2 of the two-component disk lu 2 upon the 
dimensionless wavenumber of the perturbation ka/n. Two ex- 
amples are shown: (top panel) Q g = Q s = 1.5, S 9 /S s = 0.2, 
and c/cr = 0.2, and (bottom panel) Q g = 1.25, Q s = 2.0, 
E 9 /E s = 0.65, and c/tr = 0.4. Labels near different curves mean: 
G+G - dispersion relation for two fluid components, S+G - dis- 
persion relation for stellar-fluid disk, G - for cold gas only, S - for 
stars only, and G(S) - for gas with the sound speed equal to the 
radial velocity dispersion of stars. 



Indeed, in this case ui 
we get 



K for large k. Using (EOh and m\ 
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- uj 2 / (n 2 n 2 
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(31) 



because the sum in ( |3l] ) diverges as lu — *• k. So, in realistic 
gas-stellar case this ratio goes to rather than increases. It 
is illustrated in Figure ^ where we have shown the ratio of 
relative perturbed surface densities for two cases: fluid-fluid 
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Figure 2. Plot showing ratio of the relative perturbations of gas 
and stellar surface densities for two cases: when stars are treated 
as fluid (dashed line) and when they are considered as collisionless 
(solid line) . Disk parameters here are those used in Figure ^ (lower 
panel). The distinction between two approaches is obvious: as k 
grows the amplitude of density perturbations in the gas relative 
to the density perturbations in the stars increases in the former 
case and decreases in the latter. 



and fluid-stellar with the same set of parameters. This de- 
crease in the gas density perturbation with growing k leads 
to the dominant role of stellar component in the dynamics 
of the whole system for large k, which is different from the 
JS case. 

The relative contribution to the instability per unit sur- 
face density of the component considered by JS is given by 
the ratio of the corresponding terms in the LHS of (KOI) 



7 = 



l + q 2 



(32) 



1 + q 2 R 2 

and is always larger than E 9 /E s for R < 1. In the star- fluid 
case one gets from that 

V „2 



7 : 



E s (1 + q 2 R 2 ) (1 - exp(-q 2 )J (g 2 )) ' 



(33) 



Careful study shows that for y/3/2 < R < 1 there is al- 
ways a range of small q (and k) in which 7 < E 9 /E s . Of 
course, as q grows 7 becomes greater than E 9 /E s for R < 1 
(7 — > i? _2 E 9 /E s as q — > 00 ), but for q ~ 1 the relative con- 
tribution to the instability per unit surface density is greater 
for stars if v3/2 < R < 1. Even in this case though gas con- 
tributes only several per cent smaller than the stars and in 
most astrophysically interesting cases cold material still has 
a serious impact on the stability of the system. 

Each of the criteria (|2!]) and (|3(J) produces some region 
in parameter space of the disk models in which they are sta- 
ble against local axisymmetric gravitational perturbations. 
Unfortunately, it does not seem trivial to construct an ef- 
fective analytical way of defining some effective value Q e f j 
as a function of Q g , Q a and other disk parameters so that 
the system is stable when Q e ff > 1, as it was done for the 
fluid- fluid case by Elmegreen (1995). Instead we follow the 



unstable 




Figure 3. Plot showing stable and unstable models of the two- 
fluid disk. Thick solid lines represent the marginal stability curves 
according to the star-fluid instability criterion (^) while thin dot- 
ted lines correspond to the Jog & Solomon (1984a) fluid-fluid cri- 
terion ([?o|). Each curve is labeled with corresponding ratio of the 
gas sound speed to the stellar velocity dispersion in r-direction 
R = c/cr. The part of the parameter space bounded by each 
marginal stability curve is stable for a corresponding R. One can 
see that the difference between two approaches - star-fluid and 
fluid-fluid is most pronounced at 1/Q S near 1 and I/Qg near 
where the kinetic and fluid approaches differ the most due to the 
weak influence of the gas component on the overall disk dynamics. 



approach of Jog (1996) and study the problem seminumer- 
ically. We parametrize our models here by 1 /Q s , 1 /Q g and 
R. In Figure ^ we compare the stable regions produced by 
each of the instability criteria. 

In the JS case the marginal stability curves given by 
conditions duj 2 /dk = and u> 2 = are symmetric with 
respect to the line Q a = Q g . Indeed, if these conditions are 
fulfilled and an inequality in ( ^o|) changes to equality for 
Qs — Qi,Qg = Q2 and q = q cr it, then one can easily check 
that this is also true for Q s = Q2, Q g = Qi and q — q cr jt I R 
provided that the instability condition is given by (p0|). Of 
course, this is not the case for the star-gas disk because any 
such symmetry is absent in the relation (^j|). 

As we increase R from to 1 the region of the parameter 
space occupied by stable models shrinks until R = 1. The 
further increase of R beyond 1 causes reexpansion of the 
region occupied by stable models. In fact in the fluid-fluid 
case the marginal curve corresponding to some particular 
R = Ro coincides with the curve corresponding to R = 1/Ro 
which can be directly checked using (pOj). But these models 
are likely to be uncommon since they have c > which 
seems to be unusual in real galaxies. 

The stable region in the star-gas case is in general 
smaller than that in the gas-gas case. The difference is es- 
pecially noticeable at 1/Q S ~ 1 and low 1/Q g , when the gas 
influence on the dynamics of the system is smallest and the 
stability condition is close to the one-fluid stellar stability 
criterion which is somewhat different from the fluid case. 
Nevertheless, the difference between two cases is quite small 
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for most of the parameters (especially for small Q a ) and one 
can usually use the JS stability criterion in these cases. 



5 APPLICATION TO THE SOLAR 
NEIGHBORHOOD 

In this section we apply the results derived in §^ to the 
neighborhood of the Sun in the Galaxy. 

In doing this one should realize that all the stars and 
gas in the Galaxy cannot be simply put into two distinct 
groups with some well defined velocity dispersion for stars 
and sound speed for gas. The reason for that is that the 
stars of different ages have different velocity dispersions - 
the older the star the larger its random motion. This ran- 
dom heating of stellar population is produced by the scat- 
tering of the stars by the giant molecular clouds (Spitzer 

6 Schwarzschild 1951, 1953) and/or transient spiral density 
waves in the Galaxy (Barbanis & Woltjer 1967; Carlberg & 
Sellwood 1985). 

Gas in its turn intrinsically has a multicomponent na- 
ture caused by the constant energy input from supernovae 
explosions and various cooling and heating processes deter- 
mining its thermal equilibrium (McKee & Ostriker 1977; 
Kulkarni & Heiles 1987). Recent studies have enabled us 
to distinguish 5 phases of the ISM: molecular gas in the 
form of the clouds (Scoville & Sanders 1987), cold neutral 
medium (CM) also in the form of clouds, and 3 more or less 
uniformly distributed gaseous components: warm neutral 
medium ( WNM) , warm ionized medium (WIM) (Kulkarni & 
Heiles 1987), and hot component (Savage 1987). They have 
different sound speeds and surface densities which makes 
one treat them separately. 

Unperturbed gravitational field of our Galaxy is not 
axisymmetric and it limits our analysis to some extent. It 
was also shown that the stellar distribution function can 
not always be represented by formula (Q) because nonax- 
isymmetries of the Galactic gravitational field cause ver- 
tex deviations of the velocity ellipsoid (Dehnen & Binney 
1998). Young stars of O and B types are likely not to have 
Schwarzschild distribution because they have had no time 
to be sufficiently scattered by the giant molecular clouds or 
transient spiral arms, but they probably are not important 
mass contributors in the local part of the Galaxy. The lo- 
cal approximation itself maybe questionable because most 
unstable waves have A of the order of several kpc. 

Other possible complications in real galactic disks in- 
volve the presence of cosmic rays and magnetic fields which 
could influence the gas dynamics (Elmegreen 1987). In this 
paper we are primarily interested in the purely gravitational 
aspects of the disk instability and for this reason we neglect 
them at all, though it makes our consideration less realistic 
when applied to the galaxies. For these reasons our analysis 
of the stability of the Solar neighborhood should be con- 
sidered only as mostly illustrative though bearing sufficient 
resemblance to reality. 

Following Holmberg & Flynn (2000) we split the galac- 
tic disk mass between 13 major parts: 4 gaseous and 9 stellar. 
We neglected hot component of the gas because of its low 
number density, n ~ 0.003 cm -3 , high temperature, T ~ 10 6 
K, (Savage 1987) and, consequently, large thickness which 
makes the reduction effects very important (Romeo 1992). 



For the same reason we neglect the stellar halo compo- 
nent. Parameters of the gaseous components are taken from 
Kulkarni & Heiles (1987) and Scoville & Sanders (1987). 

Surface densities of stellar components are taken from 
Holmberg & Flynn (2000). We got the radial velocity dis- 
persions based on the recent data from Hipparcos from 
Mignard (2000). Velocity dispersions of white and brown 
dwarfs have been chosen quite arbitrarily. All the parame- 
ters assumed for the mass constituents are listed in Table 
1. 
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In Figure ^ we show the dependence of u 2 upon the 
inverse radial wavelength k/2ir. Epicyclic frequency k — 36 
km s _1 kpc -1 is assumed throughout the calculation. Curve 
labeled 1 corresponds to the choice of parameters described 
in the Table 1. One can see that for all radial wavelength ui 2 
is positive that is the whole disk system is stable against lo- 
cal axisymmetric perturbations. To check how rigorous this 
conclusion is we varied some of the model parameters until 
the disk became unstable. 

Data about the surface densities and velocity disper- 
sions of the brown dwarfs (BDs) and white dwarfs (WDs) 
seem to be the most uncertain among all the model param- 
eters, so we tried to vary them first. Curve labeled 2 has all 
the surface densities as listed in Table 1 but the velocity dis- 
persion of WDs and BDs was lowered to oy = 19.2 km s . 
Only if this population is so cold can it make the system 
unstable with all other parameters being kept unchanged. It 
seems inevitable that real WD and BD populations must be 
sufficiently hotter because only young stars can have such a 
low velocity dispersion (Mignard 2000). 

Curve labeled 3 shows ui 2 — A -1 dependence for the case 
when total surface density of WDs and BDs was raised from 
10.6 M© pc -2 to 26 Mq pc -2 keeping the rest of model 
parameters unchanged. This leads to the neutral stability 
of the system but such a surface density seems to be too 
large despite large uncertainties and claims of some authors 
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Figure 4. Plot showing the dependence of uj 2 upon A -1 in the 
Solar neighborhood. Thick solid curve labeled 1 corresponds to 
the model with parameters given in Table 1. Other curves corre- 
spond to a different sets of parameters described in the text. The 
model of the Solar neighborhood described in Table 1 is stable 
against local axisymmetric gravitational perturbations. 

that such extreme values of surface density could be com- 
mon. For example, Festin (1998) found high mass density of 
the WD, ~ 2.6 times larger than we assume here, but his 
conclusions were based on a small sample of 7 sources only. 
Other authors (Ruiz & Takamiya 1995; Oswalt et al. 1996) 
claim values for WD mass density which are in agreement 
with what we take. Even more uncertainty is involved in 
determining the density of the BD. Different surveys quote 
values from 0.6 to 4 times what we assume in this research 
(Fuchs, Jahreiss, & Flynn 1998). Recent data (Reid et al. 
1999) based on a large enough sample imply mass density of 
BD 0.005 M Q pc" 3 which is about 60% of what we assume 
in our calculations. 

Finally, the fourth curve shows the dispersion relation 
for the disk with the lowered sound speed in some of the 
gas components: in CM we set c s = 5.0 km s" 1 , in WNM 
c s = 7.5 km s^ 1 , and in WIM c s = 8.0 km s _1 , with all other 
parameters unchanged. In this case disk becomes marginally 
unstable. Important thing to notice here is that small vari- 
ations in the gas sound speed can have stronger influence 
on the disk stability than large changes in the stellar veloc- 
ity dispersion, even though the surface density of the gas 
is smaller than that of stars. This is a manifestation of the 
crucial importance of cold material for the stability of the 
whole disk, which was first noted by Jog & Solomon (1984a). 
It is also easy to see that variations of the gas parameters 
produce significant change of the most unstable wavelength 
compared to the variations of the stellar parameters; it is 
reduced from ~ 2 kpc to ~ 1 kpc by that small change in 
gas sound speeds. 

The bottom line is that local Galactic disk seems to be 
stable against local gravitational axisymmetric perturbation 
even when allowance for a scant knowledge of some of the 
Galactic parameters is made, which confirms results of Jog & 
Solomon (1984b) and Elmegreen (1995) for two-fluid disks. 



6 CONCLUSIONS 

Due to the continuous interaction with the transient spiral 
structure and giant molecular clouds stars diffuse in the ve- 
locity space towards higher random velocities as their age 
increases and it was confirmed observationally (Dehnen & 
Binney 1998; Mignard 2000). It raises a necessity of consid- 
ering the dynamics of the Solar neighborhood taking into 
account its complex multicomponent structure. In this pa- 
per we studied the stability of such a system against grav- 
itational axisymmetric perturbations in the tight-winding 
limit. It is possible to derive an analytic dispersion relation 
characterizing multicomponent thin differentially rotating 
disk and study its stability. 

In doing so we distinguished between two types of the 
disk constituents: stellar and gaseous. Stellar population is 
dynamically different from fluid because stars form collision- 
less system (Binney & Tremaine 1987) while gas must be 
treated as a fluid. We demonstrated that the difference in 
the results for stability produced by two approaches is small 
for multicomponent disks in many astrophysically interest- 
ing cases. Some disk models though could be sensitive to the 
choice of the stability condition and in that case one should 
use correct criterion given by the equation (p3|). 

We apply our results to the stability of the Solar neigh- 
borhood and confirm the conclusions of the previous two- 
fluid studies that local Galactic disk is stable against ax- 
isymmetric perturbations in the WKB limit, even taking 
into account uncertainties associated with determining some 
of the disk parameters. 

Keeping in mind previous two-fluid results (JS) it is 
not surprising that relatively small variations of the gaseous 
component parameters are very important for the overall 
disk stability. Indeed, Figure [| shows that decrease of the 
sound speed of the gas by about 1 km s _1 drives instability 
of the disk, while in the case of stellar component one needs 
to reduce its velocity dispersion by ~ 10 km s _1 to produce 
the same outcome. Even though the mass of the gas in the 
disk is smaller that the stellar mass its small random motion 
makes it much more susceptible to its own self-gravity than 
the hot stellar component. 

Our study of the stability of the Solar neighborhood 
neglects a lot of physics such as magnetic fields or nonax- 
isymmetry of the Galactic gravitational field and their im- 
portance remains an open question. Measurement errors as- 
sociated with determining some disk parameters also limit 
the applicability of multicomponent criterion because the 
larger the number of constituents the larger errors get ac- 
cumulated. Future interferometric missions such as GAIA 
and SIM will probably solve this problem because of their 
anticipated accuracy. Nevertheless, even with all the sim- 
plifications and the observational uncertainties involved it 
seems that the Solar neighborhood is stable against purely 
gravitational axisymmetric perturbations. 
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